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The two-photon momentum distribution of annihilating electron-positron pairs in ferromagnetic 
nickel (Ni) was determined by measuring the spin-polarized two-dimensional angular correlation 
of annihilation radiation (ACAR). The spectra were compared with theoretical results obtained 
within LDA+DMFT, a combination of the local density approximation (LDA) and the many-body 
dynamical mean-field theory (DMFT). The self-energy describing the electronic correlations in Ni 
is found to make important anisotropic contributions to the momentum distribution which are not 
present in LDA. Based on a detailed comparison of the theoretical and experimental results the 
strength of the local electronic interaction U in ferromagnetic Ni is determined as 2.0 ±0.1 eV. 


The electronic properties of many narrow-band materi¬ 
als, such as the d-shell transition-metal series and their 
compounds, cannot be explained within a one-electron 
picture, because there exist strong correlations between 
electrons in the partially filled d band Such sys¬ 

tems are therefore better described by multi-band models 
such as the Hubbard or Anderson-type lattice model. In 
these models the local Coulomb repulsion U is assumed to 
be the dominant interaction between the electrons. The 
“Hubbard” parameter U was originally introduced for 
single-band models ii and is defined as the Coulomb- 
energy required to place two electrons on the same site: 
U = E{d^+^)pE{d^-^)-2E{dJ^). ReieEid^) represents 
the total energy of a system for which n electrons fill a 
given d-shell on a given atom. In multi-band systems U 
takes the form of an interaction matrix. 

The Hubbard model and related lattice models are able 
to explain important general features of correlated elec¬ 
trons, but they cannot describe the physics of real mate¬ 
rials in any detail. Namely, for an approach to be realistic 
it must take into account the explicit electronic and lat¬ 
tice structure of the systems. Here the LDA±DMFT ap¬ 
proach has led to great progress in the understanding of 
correlated electron materials BIMil. LDA+DMFT is 
a computational scheme where the local density approxi¬ 
mation (LDA) or the generalized-gradient approximation 
(GGA) provide the material dependent input (orbitals 
and hopping parameters), while DMFT [2|,ll2j solves the 
interacting local many-body problem originating from 
the local Hubbard interaction U and Hund’s rule cou¬ 
pling J. The results can be compared with experimen¬ 
tal data obtained, for example, by photoemission spec¬ 


troscopy. In particular, this technique measures spectral 
functions, i.e., the imaginary part of a one-particle Green 
function, and thus determines correlation induced shifts 
of the spectral weight. This allows one to estimate the 
local Hubbard interaction of a material, say, Ni. In¬ 
deed, most investigations on the electronic structure of 
Ni relied on photoemission spectroscopy lillil- Braun 
et al [m demonstrated the importance of local correla¬ 
tions in Ni by exploiting the magnetic circular dichroism 
in bulk sensitive soft X-ray photoemission measurements. 
One of the dominant correlation effects observed in the 
photoemission data for Ni is the satellite peak situated 
at 6 eV below the Fermi level EH- This feature is not 
captured by LDA, but is well explained by LDA±DMFT 
[l8|. LDA±DMFT also reproduces the correct width of 
the occupied 3d bands and the exchange splitting [l^. 
A fundamental difficulty concerning the interpretation 
of photoemission data is the fact that they involve an 
interaction of photons with matter. This makes the de¬ 
termination of a parameter such as the Hubbard U quite 
difficult. Namely, it is not only necessary to describe the 
excitation process, but also the propagation of the pho¬ 
toelectrons in the material as well as the process of de¬ 
tection itself. Therefore the experimental data, and the 
values of the interaction parameters derived from them, 
will be strongly influenced by the surface of a sample. 

In this Letter we discuss an alternative experimental 
technique to determine the local Coulomb parameter U, 
involving positrons. In contrast to photoemission experi¬ 
ments positron spectroscopy experiments measure a two- 
particle Green function. Since there are no physical pro¬ 
cesses except the electron-positron annihilation, positron 
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spectroscopy does not suffer from the above-mentioned 
difficulties faced by photoelectron spectroscopy. As a 
consequence the magnitude of U deduced from positron 
spectroscopy is much less influenced by external effects 
(e.g., surfaces), which are difficult to control. Here we 
show that by combining experimental results of the spin- 
polarized two-dimensional angular correlation of annihi¬ 
lation radiation (2D-ACAR) with LDA+DMFT compu¬ 
tations it is possible to assess the strength of the elec¬ 
tronic interactions in Ni quite unambiguously. We mea- 



Figure 1. (Color Online) Schematic principle of spin polarized 
2D-ACAR: In electron-positron annihilation the singlet con¬ 
figuration is preferred for majority or minority spin electrons 
if the magnetization of the sample is parallel or anti-parallel 
to the emission direction of the positrons. 


sured spin-resolved two-dimensional projections of the 
two-photon momentum distribution (TPMD) in mag¬ 
netic fields up to 1.1 T at room temperature. The field 
was applied parallel and anti-parallel, respectively, rela¬ 
tive to the crystallographic (110) orientation of the sam¬ 
ple which coincides with the emission direction of the 
positrons (i.e., the polarization direction of the positron 
beam). The /3+ spectrum of ^^Na, which served as 
positron source, has an end-point energy of 545 keV with 
a mean energy of about 215 keV. This corresponds to 
a beam polarization of 74 % for forward emission. Both, 
the positron back reflection in the source capsule and the 
magnetic guiding field, lower the net polarization of the 
beam, which was determined as 34.4(5) % in a separate 
experiment. Although the positron undergoes multiple 
scattering during thermalization [iMil the polarization 
of the positrons remains essentially unchanged until the 
moment of annihilation. Quantum electrodynamics pre¬ 
dicts the annihilation rate in the triplet configuration to 
be significantly smaller (1/1115) than that in the singlet 
configuration [22|, [23| . When an external magnetic field is 
applied parallel (tt) or anti-parallel (ti) to the emission 
direction the positrons will therefore annihilate predom¬ 
inantly with electrons from the majority or the minor¬ 
ity spin directions, respectively (see Fig. [T]). The result 
of a magnetic 2D-ACAR measurement can be expressed 


as 




^^eyLpiPxi Py) — Py) Py) (1) 


where N^(^^^^{px^Py) is the number of coincident photon 
counts measured with the positron spin aligned paral¬ 
lel (tt) and anti-parallel (ti), respectively, to the applied 
field. For a detailed description of the experimental tech¬ 
nique we refer to Refs. ( 2342 ^. Spin-polarized 2D-ACAR 
is one of the few experimental methods that can probe 
the momentum distribution of the electrons in the bulk 
with respect to the spin direction. It was success fully 
applied to elemental Ni [27| and other materials |28l43l|. 
In Fig. [ 2 ] we present the spin-difference of the 2D-ACAR 
measurement of Ni. In each spectrum an excess of 250 
million counts was collected, and the data were corrected 
for the momentum sampling function. Before subtrac¬ 
tion, the spectra and were normalized to an 
equal amount of counts. A renormalization due to 37 
decay was omitted since the corresponding correction in 
the case of Ni [ 2 ^ is negligible compared to the statistical 
noise. The 4-fold symmetry is clearly observed in agree¬ 
ment with the study of Manuel et al It should be 
noted that the difference spectrum (see inset of Fig. [2j) 
is anisotropic, i.e., the signal is more intense along the 
F — X — F direction than along the F — L — F direction. 



Figure 2. (Color Online) 2D-ACAR difference spectra 
^Nexp{px,Py) obtained when the magnetic field is reversed, 
with the integration along the (001) direction, px = (100) and 
Py = (010). The inset illustrate the anisotropy of the differ¬ 
ence spectra between two directions in momentum space. 

To represent the difference spectra ANexp{Px^Py) in 
k-space we apply the magnetic Lock-Crisp-West (LCW) 
back-folding procedure. In practice this amounts to tak¬ 
ing the difference of the individual spin-resolved LCW- 
projections [lH since the LCW-procedure is nothing but 
a linear transformation. AccordiM to this scheme, first 
proposed by Biasini et al. [H, [sj, effects due to the 
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electron-positron interaction, known as enhancement, are 
cancelled out. This proves to be particularly advan¬ 
tageous as no complete microscopic description for the 
enhancement problem is available. Several recipes have 
been discussed in the literature l35 -l4l| and have also 
been consequently applied for Ni 421 4^ Hence, the only 
parameter in this analysis is a uniform scale factor, which 
is obtained by fitting the amplitude of theoretical spectra 
to experimental data. 

The theoretical analysis requires the knowledge of the 
two-particle electron-positron Green function, describing 
the probability amplitude for an electron and a positron 
propagating between two different space-time points. In 
view of the fact that experiments performed on Ni |43- 
show no significant enhancement effects, we factorize 
the two-particle Green function into a product of elec¬ 
tronic and positronic Green functions. Thereby Coulomb 
interaction induced correlation effects between the anni¬ 
hilating particles are neglected. At the same time cor¬ 
relations between the electrons are explicitly included in 
the DMFT. 


Electronic structure calculations were performed with 
the spin-polarized relativistic Korringa-Kohn-Rostoker 
(SPR-KKR) method [i^. For LDA computations the 
exchange-correlation potentials parametrized by Vosko, 
Wilk and Nusair were used with a lattice param¬ 
eter of 3.52 A. To include the electronic correlations, 
a charge and self-energy self-consistent LDA-fDMFT 
scheme was employed, which is based on the KKR ap¬ 
proach [i^ and where the impurity problem is solved 
with a spin-polarized T-matrix fluctuation exchange 
(SPTF) method [13, li^- This impurity solver is fully 
rotationally invariant even in the multi-orbital version 
and is reliable when the interaction U is smaller than 
the bandwidth, a condition which is fulfilled in the case 
of Ni. In this LDA+DMFT framework the electron- 
positron momentum density Pcr(p) is computed directly 
from the two-particle Green function in the momentum 
representation. The neglect of electron-positron correla¬ 
tions corresponds to the factorization of the two-particle 
Green function in real space. In the numerical imple¬ 
mentation the position-space integrals for the “auxil¬ 
iary” Green function Gaa'{Pe^Pp) obtained within LDA 
or LDA+DMFT, respectively, are performed as integrals 
over unit cells: 


the electron-positron momentum density: 

p“(p) = -^1 dEeG",,{pe,Pp,Ee,Ep). (3) 

Integration over positron energies Ep is not required, 
since only the ground state is considered, and a' = —a at 
the annihilation. The momentum carried off by the pho¬ 
tons is equal to that of the two particles up to a recipro¬ 
cal lattice vector, reflecting the fact that the annihilation 
takes place in a crystal. Hence an electron with wave 
vector k contributes to p^{p) not only at p = k (normal 
process) but also at p = k + K, with K a vector of the 
reciprocal lattice (Umklapp process). The experimental 
spin-difference spectra ANexp{Px^Py) can be compared 
with the computed difference in the integrated momen¬ 
tum densities of Eq. [3l 

^NtheoiPx^Pv) = j dp^ [p^(p) - p^p)] . (4) 

In Fig. [3]we show the measured ANexp{Px^Py) and the 
theoretical AN^^^^{px^Py) LCW-folded difference spectra 
for different values of U. It is clearly seen that with in- 



G'^„,{Pe,Pp,Ee,Ep) = ^ J d^r j d^v' 

(r) Jm G“+(r, r', E,) (2) 

^P+ "> ’ 

where a = LDA or LDA-f-DMFT, and (pe,cr), and 
(Pp,cr') are the momenta and spin of electron and 
positron, respectively. Here is computed for each 

energy point on the complex energy contour, providing 


Figure 3. (Color Online) Experimental magnetic LCW spec¬ 
trum {center) compared to theoretical spectra computed for 
different values of the local Coulomb parameter U in the range 
from 1.4 to 2.3 eV. 

creasing U a gap opens at the L-points of the Brillouin 
zone. This gap is associated with the necks in the FS 
of Ni. Apparently LDA underestimates the density at 
the X-point, while the density near the L-point is over¬ 
estimated. In the LDA-fDMFT calculation the highest 
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density is found at the X-point, similar to the experi¬ 
mental data. However, the structure connecting the X 
and L-points is less pronounced in the experimental data 
than in the LDA-hDMFT results. Obviously the local 
interaction (see e.g. 0 ) pushes the d-bands below the 
Fermi energy, whereby the X 2 hole pocket obtained in 
LDA disappears. This also greatly changes the calcu¬ 
lated anomalous Hall effect of Ni |5Q[ l5l|. 



U(eV) 

Figure 4. (Color Online) Least square fit analysis (y^) be¬ 
tween LDA-hDMFT calculations and experimental data as a 
function of the Hubbard U for the 2D data. Higher U values 
correspond to stronger electron-electron correlations. A pro¬ 
nounced minimum of is found for U = 2.0 eV. The inset 
shows the results for the ID data. (The dotted lines act as a 
guide to the eye.) 

In order to derive the value of the local Coulomb in¬ 
teraction parameter U we performed a least square fit 
analysis of the measured data with the LDA+DMFT cal¬ 
culations (Fig. [3]). The result is summarized in Fig. 01 A 
distinct minimum in y^ is found at 1/ = 2.0 eV. The 
ID data are shown in the inset of Fig. 01 Although 
the spin difference spectra along different directions in 
momentum space are not the same, the y^ curves have 
their minima in the same region around U = 2.0 eV of 
the Hubbard interaction. The loss of information in the 
twice-integrated ID data is indicated by a larger y^ value. 
Interpolating the data in Fig. 0] with higher order poly¬ 
nomials allows us to estimate the systematic error in the 
position of the absolute minimum as TO.leV. The re¬ 
sulting value U = 2.0±0.1eV supports the conclusion 
of recent measurements of the Compton scattering pro¬ 
files [H, [ 53 I where values of 2.0eV and 2.3 eV for U were 
reported. 

To conclude, we have shown that spin-polarized 2D- 
ACAR is a powerful tool to investigate the electronic 
structure of ferromagnetic systems especially when com¬ 
bined with LDA+DMFT calculations. Here the strength 
of magnetic 2D-ACAR becomes apparent: the higher in¬ 


formation content of the two-dimensional data compared 
to twice-integrated one-dimensional data from Compton 
scattering. We found strong evidence for the existence of 
electron pockets around the X point which are related to 
electronic correlation effects. A detailed comparison of 
the experimental results with electronic structure calcu¬ 
lations in the framework of LDA-f-DMFT allowed us to 
evaluate the local Coulomb interaction (“Hubbard” U) 
as 2.0±0.1eV. 
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